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ABSTRACT 

Aims. Using the deepest (370 ksec) Chandra observation of a high-redshift galaxy cluster, we perform a detailed characterization 
of the intra-cluster medium (ICM) of WARPJ1415. 1+3612 at z=1.03, particularly its core region. We also explore the connection 
between the ICM core properties and the radio/optical properties of the brightest cluster galaxy (BCG). 

Methods. We perform a spatially resolved analysis of the ICM to obtain temperature, metallicity and surface brightness profiles over 
the 8-400 kpc radial range. We measure the following cool-core diagnostics: central temperature drop, central metallicity excess, 
central cooling time, and central entropy. Using the deprojected temperature and density profiles, we accurately derive the cluster 
hydrostatic mass at different overdensities. In addition to the X-ray data, we use archival radio VLA imaging and optical GMOS 
spectroscopy of the central galaxy to investigate the feedback between the central galaxy and the ICM. 

Results. Our spatially resolved spectral analysis shows a significant temperature drop from a maximum of 8.0 keV to a projected core 
value = 4.6 ± 0.4 keV, and a remarkably high central iron abundance peak, Zff^^= 3.60^^ gjZg, measured within a radius of 12 kpc. 
We measure MsQo=M(r < /?5oo)=2.4±0.4 Mq and a corresponding gas fraction /goj=0. 10±0.02. The central cooling time is shorter 
than 0.1 Gyr and the entropy is equal to 9.9 keV cm-. We detect a strong [Oil] emission line in the optical spectra of the BCG 
with an equivalent width of -25 A, for which we derive a star formation rate within the range 2-8 Moyr"'. The VLA data reveals a 
central radio source coincident with the BCG with a luminosity Li 4g//,=2.0x10~^ W Hz"', and a faint one-sided jet-like feature with 
an extent of ~80 kpc. We do not find clear evidence for cavities associated to the radio AGN activity. 

Conclusions. Our analysis shows that WARPJ1415 has a well developed cool-core with ICM properties similar to those found in the 
local Universe. Its properties and the clear sign of feedback activity found in the central galaxy in the optical and radio bands, show 
that feedback processes are already established at z ~ 1 (a lookback time of 7.8 Gyr). In addition, the presence of a strong metallicity 
peak shows that the central regions have been promptly enriched by star formation processes in the central galaxy already at z > 1. 
Our results significantly constrain the timescale for the formation and self-stabilization of cool-cores. 

Key words. Galaxy clusters - high redshift: observations - X-rays: Galaxy clusters - individual - WARPJ1415. 1+3612 



1. Introduction process (see Peterson & Fabian 2006' and references therein). 

In the current framework, the dominant heating source that pre- 

^ Galaxy clusters are dynamical environments hosting complex vents the overcooling of the ICM is likely to be an AGN fueling a 

_C3 ^ astrophysical phenomena, that provide us with a wealth of in- massive central black hole, leading to an outburst which can heat 

formation on the inti'icate processes that shape the cosmic lai-ge- the ICM via shocks, buoyantly rising bubbles inflated by radio 

scale structure and galaxy evolution (see review s by Fabian lobes, or dissipation of sound waves (see McNamara & Nulsen 

[19941 Rosati, Borgani & Norman 2002| Voit l2005l l. 2007 for a review). 

The intracluster medium (ICM) is the dominant baryonic The so-called cool-core phenomenon is observed in different 

component of galaxy clusters, a hot plasma emitting X-ray ra- wavebands depending on the cluster component under investiga- 

diation via thermal bremsstrahlung. High-resolution X-ray ob- tion. Among them: the ICM (X-ray), the brightest cluster galaxy 

servations show that the surface brightness of the ICM of about (BCG, optical), cold molecular gas (IR), and the central AGN 

half of the local clusters is peaked in their central regions, where (radio). The onset of star formation activity in the central galaxy 

the inferred cooling time of the gas is shorter than the typical (otherwise a passive, early-type galaxy) is a well-known mani- 

dynamical time (see Fabian 1994 for a review on cooling flows festation of this phenomenon (Crawford et al. 1 19991 McNamara 

and Hudson et al. 2010 for more recent results). However, X-ray et al. 120071 Donahue et al. 2010). However, nuclear activity in 

spectroscopy always shows a gas temperature floor, indicating the central cluster galaxy is expected to play the major role in 

that some distributed source of heating must stop the cooling regulating the cool-core thermodynamics, as suggested by the 
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observed interactions between the radio jets and the ICM and 
the scaling relations among total radio power and ICM proper- 
ties observed in the X-ray band (Birzan et al. 2008). These pro- 
cesses are shown in the spectacular combined images of nearby 
clusters in the radio and X-ray bands, where cavities and rip- 
ples in the ICM are observed to correspond to non-thermal radio 
emission (Blanton et al. '2001; Birzan et al. 2004; Wise et al. 
[200 7, Sanders & Fabian 2007, Sanders, Fabian & Tavlor |2009l 
and many others). 

There is increasing evidence from radio and X-ray observa- 
tions of local clusters that the formation of bubbles due to radio 
jets associated to the central AGN may effectively satisfy the en- 
ergy balance between cooling and heating in cool-core regions 
(Blanton |2001l Ehlert et al. l201 11 1. In this respect, the radio lumi- 
nosity of the central galaxy provides an important link between 
the black hole activity and the state of the intracluster medium. 

To date, systematic studies of the cool-core phenomenon in- 
cluding X-ray, optical and radio data to understand how the feed- 
back mechanism shapes the evolution of cool-cores have been 
limited to samples of nearby clusters (e.g., Heckman et al. 1 19891 
Birzan et al. 2004 Cavagnolo et al. 2008, O'Deaet al. l2008l 
Mittal et al. 2009 Sun 2009 McDonald et al. g OlQl*, reaching 
a redshift up to < z >=0.3 (Samuele et al. '201 li These studies 
show clear correlations involving the radio luminosity of the cen- 
tral AGN, optical emission lines (e.g. Ha), excess IR/UV emis- 
sion, and the core entropy of the ICM. 

The extension of these studies to high-redshift can provide 
important clues also on the timescales of the metal enrichment 
mechanisms of the ICM, a process that encompasses differ- 
ent phases of the cool-core phenomenon. Iron, the main metal 
locked in the ICM, is produced primarily by type la super- 
novae, while iron and other heavy elements produced by both 
type la and type II SNe hosted by the cluster early-type galax- 
ies (Renzini et al. I1993I I are eventually diffused into the ICM 
via galactic winds and ram pressure striping. According to lo- 
cal studies, the excess iron mass found in cool-cores is dkectly 
linked to the brightest cluster galaxies (Bohringer et al. l2004l De 
Grandi et al. 2004). A complete understanding of the formation 
and evolution of cool-cores should encompass all these aspects. 

Recent results indicate a moderate evolution in the bulk of 
the cool-core cluster population out to z=1.3, with an apparent 
lack of very strong cool-core clusters at high redshifts (Santos 
et al. 120081 |20l0). However, to adequately compare the char- 
acteristics of distant and nearby cool-core clusters, we require 
deep, high-resolution X-ray data of distant clusters. The study 
of cool-cores at high-redshift will not only place an upper limit 
on the formation timescale of cool-cores, but will also reveal the 
strength of the interactions between the diffuse baryons and the 
brightest central galaxy (BCG) at an epoch when the most mas- 
sive clusters are still assembling. 

In this paper we provide a detailed investigation of the ICM 
of the X-ray selected cluster WARPJ1415. 1+3612 (hereafter 
WARPJ1415) at z=1.03 (Perlman et al. l2002] l. using a recently 
acquired deep Chandra observation. This cluster was chosen as 
the strongest high-z cool-core cluster based on the analysis of its 
surface brightness properties (Santos et al. 2010). Using a multi- 
wavelength dataset comprising radio VLA data, archival opti- 
cal GEMINI-GMOS spectroscopy and imaging from HST/ACS 
F775W, and SUBARU-^M/^Wme (BVRiz) we study the central 
galaxy and its interactions with the hot plasma. The analysis pre- 
sented here is the first detailed study of the feedback process in 
a cool-core at such high-redshift. 

The paper is organized as follows: in §2 we describe the X- 
ray data reduction procedures, including the detailed spatially 



resolved spectroscopic analysis of the ICM. In §3 we obtain the 
temperature and metallicity profiles. We perform the deprojec- 
tion analysis and present the mass profile in §4. We investigate 
the surface brightness properties of the ICM in §5 and we ac- 
curately measure several cool-core diagnostics, namely cs b, the 
central cooling time and central entropy. In §6 we present the ra- 
dio and optical properties of the BCG and explore the connection 
between the BCG and the ICM core properties. Our conclusions 
are summarized in §7. 

The cosmological parameters used throughout the paper are: 
Hq-1Q km/s/Mpc, Qa=0.7 and Qni=0.3. Quoted errors are at the 
1-cr level, unless otherwise stated. 

2. X-ray data reduction and spectral analysis 

2.1. The X-ray dataset 

The galaxy cluster WARPJ1415 was detected in the Wide Angle 
ROSAT Pointed Survey (Jones et al. [T998l Perlman et al. 2002ll- 
This survey provided several high-redshift galaxy clusters based 
on serendipitous detections of extended sources in targeted 
ROSAT PSPC observations. The survey covers an area of 71 
deg^ and contains a complete sample of 129 X-ray selected clus- 
ters down to a flux limit of S ~ 6.5 x lO"'"*^ erg s~^ cm^. 

A deep, 280 ksec ACIS-S observation of WARPJ1415 was 
awarded to our group in Chandra AO 12 (PI J. Santos). An ad- 
ditional 90 ksec observation with ACIS-I (taken in 2004, PI H. 
Ebeling) was used in our analysis. The Chandra exposure time 
sums up to a total 370 ksec (obsid4163, 12255, 12256, 13118, 
13119). This dataset represents the deepest Chandra observation 
of a cluster at z ^ 1, enabling the most detailed X-ray analy- 
sis of any such high redshift galaxy cluster. Previous studies of 
distant X-ray clusters with comparable depth used both shallow 
Chandra and XMM-Newton data (see Maughan et al. 2008 ), and 
therefore do not reach the angular resolution to study the inner 
regions of the ICM distribution. In the case of WARP1415, an- 
gular resolution has a paramount relevance, since it allows us 
to resolve the density, temperature and metal distribution within 
the cool-core region on a scale as small as ~10 kpc, and hence 
to shed light on the complex physics of the cool-core region at a 
look back time of 8 Gyr 

We also examined the XMM-Newton archival data on 
WARPJ1415. We find that after removing the time intervals of 
high-background, the useful exposure time is about 17 ks for 
MOS and 14 ks for PN. Given the modest exposure time and 
the much lower angular resolution of XMM-Newton, we do not 
include this dataset in our analysis. 

2.2. Data reduction 

We performed a standard data reduction starting from the 
level= 1 event files, using the CIAO 4 . 3 software package, 
with the most recent version of the Chandra Calibration 
Database at the time of writing (CALDB 4.4.5). We ran the 
acis_process_events task to update the data calibration, in- 
cluding the tgain correctioifl Since both ACIS-I and ACIS- 
S observations were taken in the VFAINT mode, we set the 
check_v£_pha key to yes to flag background events that are 
most likely associated with cosmic rays and distinguish them 
from real X-ray events. With this procedure, the ACIS parti- 
cle background can be significantly reduced compared to the 
standard grade selection. We also apphed the CTI correction to 

' See |http://cxc.harvard.edu/ciao/why/acistgam.htmI] 
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the ACIS-I observation (taken with the temperature of the Focal 
Plane equal to 153 K). This procedure allows us to recover the 
original spectral resolution partially lost because of the CTI (see 
Grant et al. 2004). The correction applies only to ACIS-I data, 
since the ACIS-S3 did not suffer from radiation damage. The 
data were filtered to include only the standard event grades 0, 2, 
3, 4 and 6. We did not remove the acis_detect_a£terglow 
correction since all our data have a SDP version larger than 
7.4.C0. Nevertheless, we checked visually for hot columns left 
after the standard reduction, finding none. We also identified the 
flickering pixels as the pixels with more than two events con- 
tiguous in time, where a single time interval was set to 3.3 s. For 
exposures taken in VFAINT mode, there are no flickering pixels 
left after filtering out bad events. We filtered time intervals with 
high background by performing a 3cr clipping of the background 
level in the 0.5-7 keV band by using the script analyze_ltcrv. 
We did not detect any high background interval in both ACIS-I 
and ACIS-S data (nominally the removed exposure time is less 
than 0.5% for ACIS-I and less than 0.2% for ACIS-S). We re- 
mark that our spectral analysis is not affected by any possible 
residual flare unnoticed in the 0.5-7 keV band, since we com- 
pute the background from source-free regions around the cluster, 
thus taking into account any possible spectral distortion of the 
background itself. Finally, all the ACIS-S obsid were merged 
together with the merge.all tool. ACIS-S and ACIS-I are kept 
separated both for imaging and for spectral analysis, as described 
in the following sections. 

2.3. Spectral analysis 

We produced soft (0.5-2 keV) and hard (2-7 keV) band im- 
ages with the full Chandra resolution (1 pixel - 0.492 arcsec). 
We computed the cluster center in the soft-band by measuring 
the position at which the signal-to-noise ratio (S/N) within cir- 
cles of different radii is maximum. Following this procedure 
we assign the X-ray center of the cluster to RA=14:15:11.08, 
DEC=-i-36:12:03.1. This position also corresponds to the peak 
of the surface brightness emission. We remark that we did not 
detect any point source which could significantly contribute to 
the core emission, by carefully inspecting the 0.5-2.0 keV and 
the 2.0-7.0 keV X-ray images. 

The spectral analysis is performed as follows: starting from 
the cluster center, we draw rings at different radii and compute 
flie S/N in the total 0.5-7 keV band for the ACIS-S image, which 
has a signal much stronger than the ACIS-I image. In order to 
achieve a good compromise between the S/N ratio in each ring 
and the total number of rings, we set a minimum S/N > 24 within 
25", decreasing at larger radii down to 14 in the outermost ra- 
dius. This ensures that we have always more than 650 net counts 
in the 0.5-7 keV band in each ring (with the exception of the 
central bin with 500 net counts), for a total of nine rings. We 
adopted the same set of circular rings in the ACIS-I image, de- 
spite the net counts in each ring are typically less then 100. In 
the last radial bin (300-400 kpc), we used only ACIS-S, since 
there is no useful signal in the ACIS-I data. 

We extracted spectra for each ring from the ACIS-S and 
ACIS-I data separately. The response matrices and the ancil- 
lary response matrices were computed for each obsid respec- 
tively with mkacisrmf and mkwarf, for the same regions where 
the spectra were extracted. A single arf and rmf files for each 
ACIS-S spectrum were finally obtained by summing over the 
4 ACIS-S obsid with addarf and addrmf, respectively. We 

^ see |http://cxc.harvard.edu/ciao/threads/acisdetectafterglow/| 



detected 7500 (6200 in ACIS-S and 1300 in ACIS-I) photons 
within 50", corresponding to 400 kpc. 

We selected the background from empty regions of the same 
CCD in which the cluster is located. This is possible since the 
signal from the cluster has a total extension of less than 1 ar- 
cmin, as opposed to the 8' size of the ACIS-I/-S chips. The back- 
ground region is scaled to the source files by the ratio of their 
geometrical areas. In principle, the background regions may par- 
tially overlap with the outer virialized regions of the clusters. 
However, the cluster emission from these regions is negligible 
with respect to the instrumental background, and does not af- 
fect our results. Our background subtraction procedure, on the 
other hand, has the advantage of providing the best estimate of 
the background for that specific observation. 

The background is extracted (both for ACIS-I and ACIS-S 
data) from 3 circular regions nearby the clusters, up to a max- 
imum distance of 3 arcmin, and with a minimum distance of 
1 .2 arcmin. The size of each background region varies between 
50 to 70 arcsec in radius, in order to ensure the background is 
sampled from a region larger by at least a factor of 10 with re- 
spect to the outermost ring. We repeated the measurement with 
other background regions selected with these criteria and found 
no significant differences in the final results. 

The background spectra were scaled by the geometrical fac- 
tor between the background and source regions before subtrac- 
tion. Our spectral analysis is expected to be robust against back- 
ground variations and vignetting effects thanks to the high S/N 
set for our circular rings (S/N of 25 except in the two outer- 
most rings, where we have S/N=20 and 15) and to the sparse 
sampling of the background. To show this, we also compute the 
scaling factors using the ratio of the exposure maps of the source 
and background regions, and repeat the spectral analysis. In all 
cases we obtain results in excellent agreement with those ob- 
tained with the simple geometrical scaling. We remark that a 
full treatment of the vignetting effects would imply modeling 
the background by treating each component separately. We do 
not apply this procedure in our paper since it would not intro- 
duce any noticeable improvement in the final results. 

The spectra of each ring were analyzed with XSPEC vl2.6 
(Arnaud et al. 1996 ) and fitted with a single-temperature mekal 
model (Kaastra il992i Liedahl et al. 1 19951 ) using the solar abun- 
dance of Asplund et al. ( l2005lFl The fits were performed over 
the energy range 0.5 - 8.0 keV. The free parameters in our spec- 
tral fits are temperature, metallicity and normalization. The local 
absorption is fixed to the Galactic neutral hydrogen column den- 
sity (Nh) measured at the cluster position (Wilms et al. 120001 ) 
and equal to Nh = 1.05 x 10^° cm"^, and the redshift is fixed 
to z - 1.03. We used Cash statistics applied to the source plus 
background, which is preferable for low S/N spectra (Nousek & 
Shue fT989l ). 

The spectra were binned only to one photon per energy bin. 
The use of Cash statistics with background subtracted spectra 
has been used in previous works, and validated with spectral 
simulations in the case of typical AGN (see Figure 2 in Tozzi 
et al. 120061 ) ■ The same results were obtained for thermal spectra 
typical of clusters. 

The global temperature and abundance values were mea- 
sured in the region of maximum S/N in the ACIS-S data, which 
corresponds to a radius Re^, - 32" centered around the clus- 
ter position. The spectrum extracted from this region allows 



^ The values normalized to Asplund et al. (2005) are a factor 1.6 
larger than the ones normalized to Anders & Grevesse (1989) which 
have been extensively used so far in the literature. 
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us to measure temperature and Fe abundance with unprece- 
dented accuracy at such an high redshift. The best fit values are 
kT = 6.82+°-;^^ keV and Z = 0.88!° jiZo. In addition, the strong 
Ka Fe line allow us to measure the redshift with an accuracy of 
less than 1%: zp, = 1.028+° °07 (see Fig. [B- The same analy- 
sis performed with the ACIS-I data provides consistent results 
with larger uncertainties. These global values are still consistent 
with the values measured with a previous CALDB version by 
Balestra et al. ( I2007I I. using the 90 ksec ACIS-I data. Our tem- 
perature measurement is in agreement, albeit with significantly 
improved accuracy, with those found by Maughan et al. ( 2006 ) in 
their analysis of the XMM-Newton data, and by Maughan et al. 
(12008) in their analysis of the ACIS-I Chandra data. However, 
we note that previous measurements of the Fe abundance with 
Chandra gave very low values, at variance with our findings. 

We will not use the average values of temperature and metal 
abundance obtained from the highest S/N region in the rest of 
this study. 
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Fig. 1. Folded global spectrum of WARPJ1415 extracted from a 
circular region with radius of 32" in the ACIS-S data (crosses) 
fitted with a mekal model (continuous line). The redshifted Ka 
line complex of Hydrogen-like and Helium like Iron (Fe XXV, 
XXVI) is clearly visible at ~ 3.5 keV. 



3. Projected temperature and metallicity profiles 

Nearby cool-core clusters are characterized by several particu- 
lar features in the central ICM, most of which require spatially 
resolved X-ray data and are therefore challenging to measure in 
the typical low S/N observations of high-redshift clusters. In the 
following subsections, we discuss in detail the spatially resolved 
properties of the ICM in WARPJ1415. 

3.1. Projected temperature profile and map 

The most distinctive sign of a cool-core is a temperature drop 
towards the cluster center, usually down to a third/half of the 
global cluster temperature (see, e.g., Peterson et al. 120031) . A 
resolved temperature drop at high-redshift has only been ob- 
served in 3C186 at z=l.l (Siemiginowska et al. 1201 Oi l, a clus- 
ter that harbors a bright quasar in its center More recently, the 
analysis of relatively shallow Chandra data of SPT-CL J2106- 
5844, a galaxy cluster at z=l.l discovered with the South Pole 



Telescope, hints at a temperature drop from ll.O^j^ keV to 
6.5;^[ J keV in the core (Foley et al. 2011), with a drop factor 
of 1.7 ± 0.5. Similarly, the very massive merging cluster ACT- 
CL J0102-4215 at z=0.87, shows temperature variations ranging 
from 22+6 to 6.6±0.7 keV (Menanteau et al. i201 Ik coiTespond- 
ing to a drop factor of ~2, from the global, core-excluded tem- 
perature, to the cluster center. 

To measure the projected temperature profile of CL1415 we 
fitted the spectra of the rings obtained with the procedure de- 
scribed in §2. For each ring we combined the ACIS-S and ACIS- 
I data. Despite the fact that the majority of the signal is pro- 
vided by the ACIS-S data, the use of the ACIS-I spectra allows 
us to obtain slightly smaller eiTor bars with respect to the anal- 
ysis based on the ACIS-S data only. In addition, the ACIS-S 
analysis shows no significant difference with respect to the com- 
bined analysis. Therefore, from now on we refer to the combined 
(ACIS-I-i-ACIS-S) spectral analysis. 

The projected temperature profile is shown in Figure |2] left 
panel. We are able to trace the temperature profile out to 400 kpc 
with an accuracy of the order of 10%, reaching 20% in the out- 
ermost bin (1-cr error). The innermost bin has a size of ~ 20 kpc, 
while the outermost bin has a size of ~ 100 kpc. We measure a 
significant central temperature drop with Tc - 4.6 ±0.4 at < 20 
kpc, which is nearly half of the maximum cluster temperature, 
/tr = 8.0 + 1 .2 keV, reached at r = 80 kpc. The ICM temperature 
drop, from the maximum value to the core value, is measured 
to be 1.7 + 0.3. Clearly, the observed drop depends also on the 
capability of measuring the temperature in the very center of the 
cool core, which in turn depends on the angular resolution and 
on the S/N. Therefore, it is dangerous to use solely the measured 
temperature drop as an indication of the cool core strength. 

To further investigate the temperature structure of 
WARP1415, we produced a temperature map. We map a 
square region of 40" by side centered on the cluster, extracting 
circular regions spaced by 2". Each circular region has a radius 
ranging from 3" to 10" going from the cluster center to the 
outer regions. Clearly, the spectra assigned to each pixel are not 
independent to each other, so the temperature map is actually 
a smoothed map, with a smoothing length increasing with 
the distance form the cluster center. The temperature map of 
WARP1415 is shown in Figure [3l where we masked out the 
pixels where the 1 sigma error on the temperature is larger than 
50%. In addition to the cool-core, which appears smooth and 
round within a radius of 3" (the minimum smoothing length of 
the map) we notice some anisotropy in the temperature profile, 
with a difference of about -2 keV in one sector at a distance of 
about 80 kpc (the distance where the temperature is maximum). 
We will discuss further this feature in Section 6. 

3.2. Projected Iron abundance profile 

A prominent peak in the Iron distribution is always associated 
with the ICM of local cool-core clusters. The origin of this iron 
excess with respect to the almost constant value measured in the 
outer regions is ascribed mostly to Type la supernovae (e.g. De 
Grandi et al. 2004). There is no consensus yet in the literature 
on the origin of the iron peak and timescale of its buildup. While 
some studies of local cool-cores favor long enrichment times (>5 
Gyr), and describe the metal excess as a long lived phenomenon 
(Bohringer et al. 120041 ). the lack of high-z data has prevented a 
more accurate assessment of this important aspect. 

In this section, we trace the spatial distribution of iron in 
the cluster ICM out to r=400 kpc, as obtained by the spatially 
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Fig. 2. Left panel: Observed projected temperature profile of WARP1415 (red points) and the corresponding fit (red dashed line) 
compared to the deprojected measured temperature values (blue solid line). Right panel: same as left panel for the Iron abundance 
profile. The 68% uncertainty in the fitted curves is very close to the typical error bars of the measured projected quantities at similar 
radii, and it is not shown for the sake of clarity. 
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Fig. 3. Temperature map obtained in a square region 40"x40" 
centered on the cluster The colorbar indicates the temperature 
scale in keV. Pixels where the 1 sigma error is larger than 50% 
are masked in black. See text for details. 



3.3. Other metals 

In addition to Iron, alpha elements (Si, Ni, S, Mg) produced by 
core-collapse supernovae (SNII) contribute significantly to the 
enrichment of the ICM (e.g. De Grandi & Molendi 2009). A spa- 
tially resolved analysis of metals other than Iron is not feasible. 
Nevertheless, we searched for emission line of other elements 
in the spectrum extracted from the inner 32". Metal lines were 
visually explored changing manually the abundance of each el- 
ement. A given element was then removed whenever its pres- 
ence was irrelevant to the best fit spectrum. Therefore, all metals 
were initially unconstrained, and progressively we froze to 0.3 
the metals which do not contribute significantly to the folded 
spectrum. This procedure can be used to identify possible low 
S/N emission lines, by setting the corresponding element to zero 
the C-statistics increases by AC > 3. With this method, we con- 
firmed the detection of the elements Si, S and Ni at 2 cr level. The 
nominal best fit values are: Zsi = 1 .9 + 0.7 Z©; Zs = 1 .3 + 0.7 Z© 
and Z^i = 3.6 + 1.7 Z©. This is the first detection of these ele- 
ments in the spectrum of a cluster at z ^ 1 , but the lack of spec- 
tral resolution and S/N prevent us from computing a meaningful 
a/Fe profile and to derive robust constraints on the enrichment 
sources of the ICM. 



resolved, combined spectral analysis. In Fig. |2] we present the 
ICM Iron abundance profile in solar units, where the solar abun- 
dance is set to the value of Asplund et al. (12005b . We detect 
an unusually high central Fe value in the ICM of WARPJ1415, 
Zff c/Z0 = 3.60;^g g5. Even taking into account the associated 
large error bar, such a high abundance level has only been re- 
ported in the local cluster Centaurus (Graham et al. I2006I I. The 
implication of this high Iron concentration and the associated 
Iron mass are discussed in Section 6.5. We investigated also the 
Iron abundance map obtained along with the temperature map, 
but we do not find any significant anisotropy in the Iron abun- 
dance distribution. 



4. Deprojection and mass profile 

The next step in our analysis is to compute the deprojected tem- 
perature and metal abundance profiles, and then to compute total 
and gas masses. The direct deprojection of the observed profile 
with the pro j ct model within XSPEC is not feasible due to the 
errors on the temperature and the coarse binning. For the best 
exploitation of our data we proceeded as follows. First we fit- 
ted the projected temperature and Iron abundance profiles. For 
the temperature we use the functional form of Vikhlinin et al. 
(I2006l l: 



kT(r) = kTi 



(1) 



(r/w)''-' + 1 (1 + (r/nYY'b ' 
where we set b = c/0.45 (as suggested in Maughan et al. (I2007I I 
and therefore we are left with 7 free parameters. 
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This functional form is usually taken to represent the 3D 
temperature distribution, and projected to be fit to the observed 
projected temperature profile. Instead, we treat r in Eq. 1 as a 
projected radius and fit the model to the projected temperature 
profile directly. The best fitting model is then deprojected to ob- 
tain the 3D temperature profile. 

As for the Iron abundance, we adopt a double beta model 
(Cavaliere & Fusco-Feminano [T976l ): 



Z/Z0 — Z/„/Zo 



1 



(1 + (r/ri)2)/^' 



+ Z„„(/Zq 



1 



(2) 



which has 6 free parameters. Therefore we can estimate the pro- 
jected temperature and Fe abundance of WARPJ1415 at any 
given radius using equations 1 and 2 with the best-fit parame- 
ters plugged in. The next step is to deproject directly the an- 
alytical renditions of the temperature and abundance profile at 
the same time. Clearly, when deprojecting the best fit tempera- 
ture and metallicity profiles, we want to preserve the full infor- 
mation coming from the surface brightness in order to directly 
obtain also the electron volume density as a function of the 
radius. We adopt a finer binning of the surface brightness pro- 
file, with the requirement of keeping the error on the normal- 
ization of each spectrum at the level of 5-10%. This is obtained 
by selecting rings with a minimum S/N of 15. We used 20 rings 
with a width ranging from 2"to 15"for increasing radii. We as- 
signed a temperature and Iron abundance to each ring according 
to equations 1 and 2, respectively, computed at the central radius 
of each ring. We normalized each spectrum in order to have the 
same predicted net count rate as observed in the real image in the 
0.5-2 keV band. We then simulated a spectrum for each ring with 
a very high S/N. The deprojection was obtained with pro j ct, so 
that we include at the same time the effects of temperature and 
Fe abundance. In Figure 2, we show the projected data with the 
best-fit model (dashed line) and the corresponding deprojected 
profile (solid line). 

The best fit normalization of the deprojected spectra is linked 
to the electron density through the relation 



Norm — 



10- 



4n{Da{l+z)y 



nndV 



(3) 



where is the angular diameter distance to the source (cm), and 
rie and rin (cm^^) are the electron and hydrogen densities, respec- 
tively. We fitted the deprojected profile with a double beta- 
model. The electron density profile and the double beta model 
best fit are shown in Figure H] The deprojected profiles of the 
temperature and of the electron density will be used to compute 
the total hydrostatic mass, while the deprojected Iron abundance 
profile will be used to measure the Iron mass (see Section 6.5). 
In this section we only compute the total dynamical mass up to 
~ 400 kpc. 

To measure the cluster mass we adopt the usual spherical 
symmetry and hydrostatic equilibrium assumption, which leads 
to the simple equation (Sarazin[l988]l: 



M(r) = -4.0 X lO'^M0r(keV)r(Mpc) 



dlog He dlog T 
dlog r dlog r 



(4) 



The logarithmic derivative of the gas density and tempera- 
ture profile were computed numerically from our best-fit pro- 
files. When we compute the total dynamical mass profile, we 
also compute the average density contrast with respect to the 
critical density at z = 1.03, so that we can solve the equation 
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Fig. 4. Deprojected electron density profile (red points) and dou- 
ble beta model best fit show in solid blue line. The two individual 
/3 model components are shown in dashed lines. 

Table 1. Total cluster mass and gas fraction measured at three 
radii corresponding to the overdensities of 2500, 500 and 200 
with respect to the critical density at redshift z - 1 .03. 



2500 
500 
200 



(kpc) 

635+*^ 

926+^' 



3.0+°;^ X 10'-* 



fsa 



0.084+0.016 
0.10+0.02 
0.20±0.03 



M^{ris) - A4/37rr^P(.(zc/) to measure the radius where the av- 
erage density level is A. Typically, mass measurements are re- 
ported for A = 2500, 500, 200. Our results are summarized in 
Table [T] The total mass profile is shown in Figure |5] (solid line) 
with the corresponding 68% uncertainty levels. The 1 cr confi- 
dence intervals on the mass are computed by assuming that the 
relative errors on the deprojected temperature profile are equal to 
those on the projected values. Similarly, the relative error on the 
electron density is equal to half the error on the spectral normal- 
ization (therefore at the level of 5% as described above). This 
is a consequence of the fact that we do not deproject directly 
the observed profiles, but rather the analytical fitting formulae, 
therefore we do not introduce additional noise. Clearly the pro- 
cedure is correct as far as the fitting formulae accurately repro- 
duce the projected profiles. This is a fair assumption, given the 
large number of free parameters, which allow us to neglect sys- 
tematic errors associated to the assumed analytical models. 

We remark that we measured temperature and density up 
to 400 kpc, which is the upper bound of the last bin with sig- 
nificant signal. Beyond this radius, all the derived quantities 
are obtained by extrapolating the best-fit profiles. As an ex- 
ample, if we assume kT - const for R > 400 kpc, instead 
of extrapolating the analytical temperature profiles, we mea- 
sure M2(K) = 3.7+°-^ X IO'^'Mq (a value 25% higher). The weak 
lensing analysis of WARPSJ1415 has been recently published 
in Jee et al. ( 1201 11 1 using HST/ACS data, yielding a total mass 
M,o,(r < 1.09 Mpc) = 4.7+J-4 x IO'^^Mq, that is consistent with 
the X-ray hydrostatic M2(K) for kT = const above 400 kpc, and 
marginally consistent with that obtained with a straight extrapo- 
lation. We defer to a forthcoming paper the detailed discussion 
of the X-ray mass profile of WARPJ1415, as well as independent 
measurements of the total mass at different radii from a dynami- 
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cal analysis of all cluster members and the modelling of a strong 
lensing system. 



1 arcsec 1 00 



lO" 



10'3 













1 1 . ■ .. ' 















100 1000 
r kpc 

Fig. 5. Total mass for WARPJ1415 from hydrodynamical equi- 
librium (red continuous line) and gas mass (green dashed line). 
The two dashed lines above and below the total mass profile 
show the 1 cr confidence interval. The profiles are shown with 
thick lines up to 400 kpc. Thin lines at r > 400 kpc show the ex- 
trapolation beyond the X-ray detectable emission. Dashed verti- 
cal lines corresponds to /?25()0, ^500 and /?2()0, from left to right. 
The red data point refers to the mass derived with the global 
cluster temperature, whereas the blue point corresponds to the 
weak lensing mass (Jee et al. 1201 11 1. 



combination of their intrinsic small angular size and X-ray data 
that usually provides only low photon statistics, make it diffi- 
cult if not impossible to distinguish between a single- or a two- 
component y6-model. 

We fitted the radial profile of WARPJ1415 using both the 
single- and double-/? model approximations (see Fig. |6l). Given 
the high-quality of our data, we can measure a significant im- 
provement in the fit using the double-beta model, with respect to 
the single-beta model, not only qualitatively but also statistically 
- see Table 2 for a list of the fit parameters. 



5.2. Surface brightness concentration, csb 

In Santos et al. (I2008I I. we defined the phenomenological param- 
eter cs B that quantifies the excess emission in a cluster core by 
measuring the ratio of the surface brightness (SB) within a ra- 
dius of 40 kpc with respect to the average SB within a radius of 
400 kpc: Csb ^ SB{r < 40kpc)/SB{r < AOQkpc). This simple 
parameter is particularly useful when dealing with low S/N data. 
The inner aperture with a physical size of 80 kpc corresponds to 
the typical size of cool-core clusters regardless of their redshift. 
We stress that we use a physical radius instead of a scaled radius, 
since the cool-core phenomenon is non-gravitational in nature, 
and therefore the self-similar scaling relations are not appropri- 
ate in this case. 

In Santos et al. ( 1201 Oi l, we measured csb=0. 144+0. 016using 
the 90 ksec ACIS-I observation. The new, more accurate value 
measured in these deep observations, c^g^O. 150+0.007, is con- 
sistent with the previous one. 



We also computed the total baryonic mass contributed by the 
ICM. The mass in baryons is obtained by directly integrating the 
electron density shown in Figure|4l assuming the rtH^nJl .2. The 
measured gas fraction at rsoo is 0.10 + 0.02 (see Table 1 for fh 
measured at different radii), typical of other distant clusters (see 
e.g. Ettoriet al. l2009] l. 

5. Surface Brightness properties, cooling time and 
entropy 

The most immediate observational signature of the presence of 
a cool-core is a central spike in the surface brightness profile of 
a cluster In this section, we use the cluster surface brightness 
properties to derive several cool-core diagnostics. 

5.1. Surface Brigiitness profile 

We computed the azymuthally averaged surface brightness pro- 
file out to r-l Mpc using the vignetted corrected ACIS-S merged 
image. The isothermal jS-model proposed by Cavaliere & Fusco- 
Femiano (1976) is often used as a simple description of the X- 
ray surface brightness profiles of galaxy clusters. 

5(r) = 5o(l+(r/r,)V^"' + C (5) 

where S[), rc, fi and C are the central surface brightness, core 
radius, slope and constant background, respectively. The fitting 
procedure is based on a Levenberg-Marquardt least-squares min- 
imization and constrains the parameter fi to the range 0.4< fi 
<1.0. Local cool-core clusters require a double j6-model to cap- 
ture the central emission peak. For distant clusters however, the 



5.3. Core luminosity excess 

Studies of nearby galaxy clusters have estimated that the ICM 
core contributes to about 25% of the total cluster luminosity 
(Peres et al. 1998 , Best et al. ,200 7). The method used to compute 
this luminosity boost was to measure the ratio L{r < rc„oi)/Lhoi. 
We argue that this is not the most efficient way to isolate the con- 
tribution of the cool-core to the luminosity because a fraction of 
the flux in L(r < rcooi) comes from the bulk of the cluster. 

To quantify the excess luminosity due to the cool-core, we 
simply computed the ratio of the flux enclosed by each of the yS- 
model components (core and outer part) at r=40 kpc (see Table 
1). This is the radius used in csb and is also very close to the 
crossing of the two j6-model components. Following this ap- 
proach, we estimate that the core contributes 4 times more flux 
relative to the bulk of the cluster luminosity, at the cooling ra- 
dius. 

The cluster soft-band luminosity within r=400 kpc is equal 
to (2.83+0. 14)xl0'*"^ erg s"'. We extrapolate this observed lumi- 
nosity to a radius of 1 Mpc using the single fi model, and obtain 
L(r < 1 Mpc)=4.0xl0'*'^ erg s"' (note that 1 Mpc is about Rom 
according to our Table 1). 

An accurate assessment of the excess core-luminosity has 
implications for the completeness of X-ray selected cluster sam- 
ples. In principle, in a purely flux-limited sample, cool-core clus- 
ters are preferentially selected with respect to non cool-core 
clusters with the same mass, thanks to their higher Lx. A proper 
treatment of this effect will be relevant to properly measure the 
evolution of cool-cores in future X-ray surveys (see Santos et al. 
EOTTT l. 
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Table 2. Single- and double-yS model fit parameters: (1) central SB in cts/s/arcmin^; (2) slope /3l; (3) core radius 1 in kpc; (4) 
reduced ;t'2; (5) (6) and (7) correspond to the second model component: central SB; slope /3 2 and core radius 2, respectively. Errors 
are not presented when the parameter value is at the limit imposed by the fitting procedure. 



Fit 


SOI 






rci 


X2 


S02 


/32 


rc2 




(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


l-fi 


0.135 


tO.007 


0.51 


50±2 


12.27 








2-/3 


0.457 


tO.067 


0.60+0.15 


14±4 


3.89 


0.044±0.011 


0.74±0.08 


150±27 




Fig. 6. Radial surface brightness profile of WAPSJ1415 (black dots). (Left) Single-beta model fit and (Right) double-beta model fit 
in red solid line. The two components of the double-yS model are shown in blue dotted lines. 



5.4. Entropy 

Clusters are usually divided in cool-core and non cool-core 
on the basis of the central value of their entropy and cooling 
time. However, there are no physical reasons to expect the ex- 
istence of two distinct populations, in fact the data rather show 
a transition between these two categories. The specific entropy, 
K(r) - kTnf'^^, is a widely used quantity to describe the thermo- 
dynamical history of the ICM. Previous studies at low-redshift 
have shown that cool-core clusters have a significantly lower 
central entropy (Kc < 30 kev cm^) than non cool-core clusters 
(Cavagnolo et al. 2008) . Using the deprojected temperature and 
gas profiles as described in Section 4, we obtain the entropy pro- 
file in Fig.|2] 

In order to compare the central gas entropy of WARPJ1415 
with the quoted values at low redshift, we have to keep in mind 
that, since there is no defined radius to perform this measure- 
ment, the local measurements are done in very small radii (< 1 
kpc), whereas our data do not allow us to go below r=8 kpc. 
Thus, the central entropy measured in the innermost bin with 
r=8 kpc, is ^"^=9. 9+2.0 keV cm^, which immediately places 
WARPJ1415 in the cool-core regime. 



1000 




10 100 
kpc 



Fig. 7. Entropy profile of WARPJ1415 computed from the nu- 
merically deprojected temperature and density profiles. 



Note that the central entropy value obtained from the depro- 
jected profiles refers to a bin centered at 4 kpc, and therefore it 
is nominally computed at a scale below the actual resolution of 
our data. We also computed the central entropy at a radius r - \2 
kpc as in the first spectroscopic bin, where we actually measured 
the temperature in the projected data (see Fig.|2l left panel). At 
r = 12 kpc the projected temperature is kT - 3.8 ± 0.3, and this 
gives a value of Kc(r - \2kpc) - 20.9 + 2.7 kev cm^. This value 
represents our conservative estimate of the central entropy and it 
is consistent with a linear relation K(r) oc r for r < 40 kpc down 
to r < 10 kpc. 



5.5. Cooling time 

The ICM central cooling time is the quantity most often used to 
characterize and quantify the cool-core strength of a cluster (e.g. 
Hudson et al. 1201 Oi l. Without a heating source to compensate 
radiative cooling, the ICM will radiate its thermal and gravita- 
tional energy on a timescale tcoi - pIKj - ^)nenHMT)] < 1 
Gyr (Fabian & Nulsen [T977l l, where p is the gas pressure, A(r) 
is the cooling function and y is the ratio of specific heats of the 
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Table 3. X-ray cool-core estimators. The quoted temperature 
and metallicity values correspond to the projected quantities, 
whereas the central cooling time and central entropy values are 
based on deprojected quantities. 



Cool-core estimator 




T _ T _ T 

core -* max -* av§ 

CSB 

^cool,c 

K, 


4.6 - 8.0 - 5.7 [keV] 
3 60+' 
0.150±0.007 
0.06+0.01 [Gyr] 
9.9±2.0 [keV cm^] 



gas. Adopting an isobaric cooling model for the central gas, tcooi 
can be computed as; 

2.5n^T 

tcool(r) = ...rr ^ (6) 

where MT, Z), Ug, and T are the cooling function, gas number 
density, electron number density and temperature respectively, 
with Yig-XSrie (Peterson & Fabian 2006). Local cool-core clus- 
ters are defined by a central cooling time much lower than the 
Hubble time (typically tcooi,c ^1 Gyr). 

An accurate measurement of the ICM cooling time profile 
of distant galaxy clusters has been unattainable up to now, due 
to the use of a single ICM temperature and the average electron 
density instead of the resolved profiles. With the current data, we 
are able to trace the precise cooling time profile of WARPJ1415, 
using the resolved temperature, density and metal abundance 
profiles. The cooling function A{T,Z), was computed using the 
cooling tables from Sutherland & Dopita (I1993I I that account for 
a varying metallicity. 

By proceeding similarly to the computation of the central 
entropy, we obtain a central cooling time tcooi ranging from 
0.06+0.01 Gyr at r = 8 kpc to 0.23 Gyr at r = 12 kpc (the central 
bin of our spectroscopic analysis). This accurate measurement is 
a factor 14 lower than the upper limit of tcooiir <20 kpc)=3.4 Gyr 
obtained by Santos et al. (1201 Oi l, using the archival data and the 
global cluster properties. This result clearly demonstrates that, 
unless deep, high-resolution data of high-z clusters are available, 
the central cooling time can only be taken as an upper limit. 

In Table 2, we summarize the X-ray cool-core properties of 
WARPJ1415 using the various cool-core diagnostics described 
above. 

6. Connection between ICM core properties and the 
BCG 

In this paper, we also explore the connection between the ra- 
dio and optical properties of the brightest cluster galaxy with the 
ICM cool core properties. Brightest cluster galaxies hold a spe- 
cial place in the history of galaxy formation and evolution. They 
are the most massive galaxies in the Universe, and are thought 
to have developed through mergers as expected in a hierarchical 
assembly model. They are located at the bottom of the potential 
well of massive clusters and are often coincident with the X-ray 
peak emission of the hot gas permeating galaxy clusters. The 
role of the BCG in shaping the ICM properties in the core as a 
function of cosmic epoch remains unclear, but it is expected to 
be an important factor in the ICM evolution, playing an impor- 
tant role in the thermodynamical equilibrium of the cool-core. 

In the center of cool-core clusters, star formation arises as a 
result of the cooling process of the ICM (Crawford et al. ll999l l. 



+36-12'30 



— +36'12' 



+35'11'30 




14hl5ml5s 14hlBml0s MhlBmOSs 

R.A. (J2000) 



Fig. 8. Subaru-Suprime BVR color image of WARPJ1415 cover- 
ing an area with 2'x2' (~ 1 x 1 Mpc), centered at RA=14:15:11, 
DEC=+36: 12:03. North is up and East is to the left. X-ray con- 
tours (with levels [3,5,10,20,30,50] cr above the background and 
smoothed with a gaussian kernel with FWHM=5") are overlaid 
in green. The 15"xl5" inset shows the core with a zoom factor 
of 2. 



with typical star formation rates in their BCGs on the order of 
1-10 Mo/yr (O'Dea et al. '2008). Interestingly, this connection 
may not be strictly localized in the cluster center, as indicated by 
McDonald et al. (1201 Oi l, who found a correspondence between 
the spatial location of Ha emission regions with the X-ray mor- 
phology (un/disturbed) of the ICM core region in nearby sys- 
tems. 

The radio emission of the central galaxy is related to the ac- 
cretion process onto a supermassive black hole. The incidence 
rate of radio sources in the center of cool-core clusters is at least 
70% (Best et al. 2006, Mittal et al. 2009; Dunn et al. 20081, 
whereas fewer than 30% of non-CCs host central radio sources. 
Several studies have established a link between the radio lumi- 
nosity of these sources with the X-ray and optical properties 
of their host clusters. In addition, other morphological features 
such as radio lobes or jets are occasionally detected. At redshift 
greater than 0.5, interactions between the central galaxy and the 
ICM remain largely unexplored. 

The BCG of WARPJ 1 4 1 5 is very large, luminous and unusu- 
ally massive, with a reported stellar mass of 2x10'^ Mo (Fritz 
et al . 120091 1. The color image obtained with SUBARU-^Mpnme 
Cam BRZ (Fig. |8]l shows that the cluster central galaxy stands 
out as a large, massive red galaxy, and we confirm that the BCG, 
the central radio source and the peak of the X-ray emission, are 
spatially coincident well within 1". In this section, we analyze 
the radio and optical properties of the central brightest galaxy, in 
order to investigate the feedback mechanism between the BCG 
and the ICM. 

6.1. Central radio galaxy and extended structure 

WARPJ1415 was observed with the Very Large Array (VLA) of 
the National Radio Astronomy Observatory in 2002 and 2003 in 
the A and B configurations at 1.4 GHz (proposal code AP439). 
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Fig. 9. (Left) Chandra full-band image with VLA radio contours overlaid in green. (Middle) Residual Chandra soft-band image 
after subtraction of the best-fit single yS-model. A lack of X-ray emission is observed in sector A with respect to the average surface 
brightness. This asymmetry is quantified by comparing the average emission in the circular regions labeled 1 and 2 (see text). (Right) 
VLA image of the radio source coincident with the cluster core. A pronounced extended feature with an extent of 80 kpc is seen in 
the NW direction. All images have a size of I'xl'. The Chandra images are shown in logarithmic scale with a Gaussian smoothing 
of 2". 



The combined radio map, with one hour exposure, has a resolu- 
tion of 2" and a noise level of 0.016 mJy/beam. It shows a bright 
radio source coincident with the X-ray centroid of the cluster and 
the BCG. Thanks to our high resolution the radio source is re- 
solved in a bright nuclear emission (3.59+0.02 mJy) and a fainter 
one-sided structure in PA ~ 70° with an extent of ~ 10", corre- 
sponding to a physical size of 80 kpc at the cluster redshift (see 
Fig-0- The flux density of this structure is (0.71+0.05) mJy and 
therefore the total flux density of the source is 4.3 mJy. Using 
low resolution NVSS data we estimate a total flux density of 
about 5 mJy, suggesting that in our high resolution image all the 
radio source is visible. 

The source radio power is in the range between FanarofF- 
Riley FR 1 and FR 11 radio galaxies, with a total radio power 
of 2.0x10^^ W Hz-i (vLv=2.8xl04i erg s"'). According to the 
low-redshift study of Sun et al. ( i2009) . aU BCGs with a radio 
AGN more luminous than 2x10^^ W Hz ' at 1.4GHz are found 
to have X-ray cool-cores. We also searched for the location of 
WARPJ1415 in the BCG radio power vs. Kc correlation for clus- 
ters with z < 0.2 presented in Cavagnolo et al. (2008) (Figure 2 
of that paper). We find that our cluster falls in the area under the 
threshold value <30 keV cm^, and with vLy >10'*" erg s"\ 
that characterizes nearby cool-core clusters. 

The extended, asymmetric radio structure could be inter- 
preted as a mildly relativistic one-sided jet. However this is un- 
likely, since in relatively low-power radio galaxies jets are rel- 
ativistic only on a few kpc scale, therefore with a size of ~ 80 
kpc we should expect to detect radio emission also on the other 
side of the core. No radio lobe seems to be present since the total 
radio flux density at low resolution is not too far from our high 
resolution result, and furthermore the structure is resolved also 
transversally. We may interpret this radio morphology as a tail- 
like structure due to a strong interaction of extended lobes with 
the surrounding medium. 

Assuming equipartition conditions, we can estimate that in 
the extended emitting region (~ 40 kpc) the magnetic field is ~ 
6 //Gauss and the minimum non thermal energy is ~ 3 x 10 
erg cm"-'. 



Deeper and higher angular resolution radio data are required 
to confirm the morphology of the extended radio emission and 
to provide a direct evidence of a jet-lCM interaction. With the 
present data, we can only report the tantalizing hint of the feed- 
back mechanism in action in WARPJ1415. However, the high 
radio power of the central AGN, and the moderate residual star 
formation estimated in the BCG (see next section) add further 
evidence in favor of a scenario in which the energy released by 
the radio loud AGN activity is able to stop the central cooling 
flow. 

6.2. X-ray cavitives? 

X-ray cavities originated by outflows of a central radio source 
are occasionally detected in the ICM of local clusters (e.g. 
Birzan et al. .2004. Fabian et al. 2006) . The detection of such 
bubbles is challenging even at low redshift due to their low sur- 
face brightness contrast and small sizes that range from <1 kpc 
to ~40 kpc (larger cavities are very rare). Detection rates are of 
the order of 20-25% in the Chandra archive (RafFerty et al. 2008) 
and in the flux-limited B55 cluster sample (Dunn et al. I2005I I. 
Nearby strong cool-core clusters have a much higher incidence 
rate, reaching 70%. 

Even though we have a deep, high-resolution observation, 
distant clusters have a small angular size which makes the de- 
tection of X-ray cavities very difficult. We searched for X-ray 
bubbles and/or surface brightness asymmetries that could be as- 
sociated with the extended radio emission. In order to do this 
we explored several techniques to improve the image contrast, 
namely, by subtracting both single and double j6-models to the 
original Chandra soft-band image, and by creating an unsharp- 
masked image. In the residual images it was difficult to see a 
convincing indication of a cavity. However, there is a clear asym- 
metry in the residual image obtained by subtracting the single-/? 
model (see Sect. 5.1) to the X-ray data, with a significant lack of 
X-ray emission in the upper half of the SB map (labeled sector A 
in Fig.|9]i, in the direction of the extended radio emission (green 
contours). We quantified this asymmetry by measuring the num- 
ber of net photons in two circular regions with a radius equal to 
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10"in the exposure-corrected Chandra image (see circle 1 and 
2 in the middle panel of Figure 9). The two numbers should be 
consistent with each other if no asymmetry is present. Instead, 
we measure 192 + 14 photons in circle 2 and 153 + 12 photons in 
circle 1, showing that region 1 (the one in the direction of the jet) 
is 25% less luminous with respect to region 2 at a 2-cr confidence 
level. 

6.3. Equivalent width of the [Oil] emission line 

Optical emission lines (Ha at 6563 A) are a relatively dust- 
independent measure of recent star formation. In alternative, the 
[Oil] (A 3727) emission line is a good proxy for Ha. Brightest 
cluster galaxies are usually red, early-type galaxies, undergo- 
ing passive evolution. Hence, star-formation is seldom found in 
these central, massive, cluster galaxies. In spite of this, high star 
formation rates have been found in BCGs of local clusters with 
central cooling times shorter than ~ 0.5 Gyr (Rafferty et al. 120081 
Cavagnolo et al. l2008l l. most likely caused by a residual cooling 
flow. 

The brightest central galaxy of WARPJ1415 was observed 
with the Gemini Multi-Object Spectrograph (GMOS-N, Hook et 
al. l2004] i in 2003 (PI Ebeling). We reduced the archival data with 
IRAF procedures and obtained the final spectrum, with a scale 
of 1 A/per pixel and a resolution R x 1000, as the result of the 
median stacking of 10 spectra with a total exposure time of 5.3 
h and a S/N~5. We expect no major slit loss since the average 
seeing is 0.7" and the width of the slit is 1". 

A significant, broad [Oil] emission line is found at the cluster 
redshift. The broadening of optical emission lines originated by 
warm ionized gas has been often observed in the BCGs of cool- 
core clusters and is likely a Doppler effect caused by motions in 
the center of the galaxy (Heckman et al. 1989). 

We measured the equivalent width (EW) of the [Oil] line 
defined by 



EW ^ 




where Fc is the continuum flux and F,i is the flux of the emission 
line. The rest-frame EW of the [Oil] line in the central galaxy of 
WARP1415 is -25+3 A. Our result goes against the interpreta- 
tion by Samuele et al. (201 1) of a strong decrease of star forma- 
tion activity with redshift, based on the EW of the [Oil] line in 
77 BCGs selected from the 160 SD survey. In that work, BCGs 
with an EW stronger than -15 A are not present. As mentioned 
in Santos et al. (1201 Oi l, we suggest that cool-core clusters might 
be under represented in the 400 SD sample (and hence its subset, 
the 160 SD sample), which might also explain the lack of strong 
[Oil] lines in the clusters' BCGs. 

6.4. Star-formation rate of the BCG 

We estimate the intrinsic [Oil] Une luminosity assuming E(B- 
y)=0.3, using the following relation from Kewley et al. (120041) : 

L[OII]i(erg s"') = L[OII]o x lO^ "^ (8) 

where L[OII]„ is the observed [Oil] luminosity. The [Oil] line 
falls in optical i-band at the cluster redshift, therefore we used 
the high-resolution HST/ACS F775W (/'-band) archival observa- 
tions to determine the luminosity of the [Oil] line. We scaled the 
observed spectrum to the F775W magnitude of the BCG mea- 
sured within a radius of 0.75", in the interval encompassed by 
the F775W filter. 
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Fig. 10. Observed optical GMOS spectrum of the BCG of 
WARPJ1415 showing a prominent, broad [Oil] emission line. 

To derive the galaxy star formation rate, we applied the 
modified Kennicutt law corrected for reddening as presented in 
Kewley et al. (120041 1. considering a mass range of 0.1-100 M© 
for a Salpeter (Salpeteril955) initial mass function: 

SFR [Oil] (Msyr-^) = (6.58 + 1.65) x lO-'^^L[Oll]. (9) 

Depending on the assumptions made on the amount of dust, we 
provide a range of SFR. In a dust-free scenario, we estimate 
SFR=2.2 Mg/yr, instead, if we assume a reddening correction 
of E(B-V)=0.3, a value typically used in star-forming galaxies 
(e.g. Lemaux et al. 2010), we obtain SFR=8.3 M^yrK 

Alternatively, the [Oil] line emission in red-sequence galax- 
ies of high-redshift clusters (z ~0.9) may be caused by an AGN 
(Lemaux et al. 2010). Although we cannot entirely rule out emis- 
sion from an AGN contributing to the estimated star formation 
rate in WARPJ1415, the fact that we do not detect an X-ray 
point source coincident with the BCG supports the assumption 
that the [Oil] emission can be entirely ascribed to star formation 
processes associated to the residual cooling flow in the core. In 
a forthcoming paper, we will further explore the properties of 
the BCG and the cluster galaxy population, using a rich optical- 
IR dataset, including the photometry used here from HST/ACS, 
SUBARU/Suprime in addition to Spitzer-IRAC. The use of in- 
frared data will allows us to firmly disentangle a possible AGN 
contamination to the SF rate diagnostics. 

6.5. Fe excess and the ICM metal enrichment process 

The interplay between the BCGs and their host clusters were 
first studied by Edge & Stewart ( 119911 ). They found a correla- 
tion between the optical luminosity of the BCG with the X-ray 
luminosity and hot gas temperature of its host cluster This work 
was expanded by De Grandi et al. (2004), where the excess iron 
mass in the core region of the ICM has been shown to correlate 
with the optical and NIR luminosity of the BCG, suggesting a 
co-evolution between the stellar content of the galaxy and the 
ICM metal abundance. 

We measured a total Fe mass Mf, = 6.2tf-9 X 1O''M0 (within 
/?20o) using the best-fit abundance profile (see Section 4). For the 
purpose of our study, it is more interesting to measure the Fe 
mass excess, M"^, following the definition of De Grandi et al. 
(120041) . This is simply the residual iron mass within 200 kpc, ob- 
tained by subtracting the almost constant abundance value mea- 
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sured at R > 200 kpc (0.25Zo in units of Asplund et al. l2005T l. 
We measured M^*J=1.8^Qg x lO^M© which corresponds to about 
30% of the total Fe mass. This number is larger than the average 
value of 10% found in nearby cool-core clusters studied in De 
Grandi et al. (2004), but the difference is not significant due to 
the large scatter in the observed local values. In fact, we find that 
our cluster lies in good agreement with the correlation between 
the iron mass excess measured in local cool-cores and their tem- 
perature. 

To probe the connection between the optical properties of the 
central galaxy and the metal enrichment of the ICM, De Grandi 
et al. (120041 1 found a correlation between the absolute optical 
magnitude of their local BCGs computed with the Rc-band with 
the excess iron mass. Assuming that the iron excess is originated 
entirely by the BCG, this relation may imply that the efficiency 
of the metal transport mechanisms from the galaxy to the ICM 
may be approximately the same in all clusters. 

In order to check whether WARPJ1415 is consistent with 
this correlation found in local cool-cores, we measured the aper- 
ture magnitude of the BCG using the Suprime Rc-band. The k- 
correction was computed with a SWIRE elliptical template of 5 
Gyr, generated with the GRASIL code (Silva et al. [T998b . We ob- 
tain an absolute magnitude of -23.8 mag, therefore WARPJ1415 
falls in the expected Moptwai - M"^ relation for local cool-core 
clusters (see Fig. 10 of De Grandi et al. 2004). In the assump- 
tion that the excess iron mass is entirely due to Type la SNe and 
stellar mass loss in the BCG, we can immediately place an up- 
per limit to the time-scale needed to build up the iron peak. We 
can safely assume that the bulk of the stars in the BCG formed 
at about z >2 (Renzini l2006l ). We also consider z = 3 as upper 
limit to the epoch when most of the stars in BCG have formed, 
as indicated by the observation of major star-formation episodes 
in massive spheroids at z > 2 (Daddi et al. 120041) or stellar pop- 
ulation modeling of high-z BCGs (e.g. Rosati et al. 2009 ). This 
corresponds to a lookback time of 10.3 Gyr (11.5 Gyr at z=3) 
in the adopted cosmology. Since WARPJ1415 lies at a lookback 
time of nearly 7.9 Gyr, the build up of its metal peak must have 
happened on a timescale of ~2.4 Gyr with an upper limit of 3.6 
Gyr. Our results point towards times scales shorter than previ- 
ously claimed (>5 Gyr see Bohringer et al. 2004) to develop 
a supra-solar iron abundance peak. This would imply a SNe la 
rate larger than expected at z > 1 in the BCG, an evidence which 
should be compared with the recent first estimate of Type la SN 
rate in high-z clusters (Barbary et al. 120 12|). In a forthcoming 
paper (Tozzi et al. in prep.) we will investigate the far reaching 
consequences these data have on chemical enrichment models, 
specifically how one can constrain the time scale of the iron pro- 
duction rate through SNe la, and the relative contributions of 
SNe la, stellar mass loss and SNII. 



7. Conclusions 

In this paper we presented a unique Chandra observation of a 
distant galaxy cluster These data enabled a spatially resolved 
analysis of the ICM, with the main goal of studying the prop- 
erties of the cluster cool-core, and trace a connection with the 
central galaxy, using optical and radio data. We measured the 
following X-ray properties: 

- a significant temperature decrease towards the center: 
rc=4. 6+0.4 keV, is 2.2 keV lower than the global tempera- 
ture, r=6.8 keV and 3.4 keV lower than the maximum tem- 
perature at r=90 kpc, r=8.0 keV; 



- a surprisingly pronounced metallicity peak, Zpc,c- 
3.6O^,jg5Z0 with a corresponding excess iron mass 
M^;"=1.8+JJ;^ X IO'^Mq; 

- a central surface brightness excess modeled by a double-/? 
model and quantified hy csb - 0.150 + 0.007. The contribu- 
tion of the cool-core to the X-ray luminosity within r = 40 
kpc is 4 times the value obtained by extrapolating the sur- 
face brightness profile from outside the cool-core towards 
the center; 

- the central cooling time, fcoo;,c=0.06 (0.23) Gyr and the cen- 
tral entropy, K,^9.9 (20.9) keV cm^ at r = 8 (12) kpc (the 
values in parentheses correspond to a more conservative es- 
timate); 

- using the measured temperature and density profiles, the 
cluster total X-ray mass, under the assumption of hydrostatic 
equilibrium, is M2()o=3.O!jJ-^xlO''*M0. 

Using VLA high-resolution data we detected a source coin- 
cident with the BCG with a radio luminosity Li.4g//j.=2.0x10^^ 
W Hz A faint, one-sided structure with an extent of 80 kpc 
is seen in the north-west direction, where a significant lack of 
X-ray emission was found. Furthermore, the analysis of optical 
spectroscopy of the central galaxy shows a broad [Oil] emis- 
sion line, with a moderate to strong equivalent width of -25 A, 
corresponding to an associated star formation rate in the range 
[2.3-8.3] Mo/yr. 

We were also able to trace a connection between the radio 
and optical properties of the central galaxy with the X-ray cool- 
core. By comparing the correlations among the radio luminosity, 
SFR and Kc measured with low-redshift clusters with our results, 
we confirm the same feedback mechanism at work in the core of 
WARPJ1415. 

We were also able to investigate the connection between the 
radio and the optical properties of the central galaxy with those 
of the ICM properties in the core. In particular, the ratio between 
radio luminosity, Kc magnitude and the estimated SFR in the 
BCG of WARPJ1415 is consistent with what is found in low- 
redshift clusters. These findings suggest that the feedback mech- 
anism at work in the core of WARPJ1415 is of the same kind 
and intensity of that observed in local clusters. 

The prominent Fe peak indicates that the metal enrichment 
mechanisms by type la supernovae (SN la) and star formation in 
the BCG happened on a short timescale (given the lookback time 
of 7.8 Gyr), and/or that the transport processes that drive away 
the metals to the outskirts (e.g. galactic winds) were not efficient 
to smear out the Fe excess. 

The analysis of the Chandra data shows that WARPJ1415 at 
z-l has all the classical features that characterize nearby cool- 
core clusters. These observations enabled the most detailed anal- 
ysis of the ICM of a z > 1 cluster and our results highlight the im- 
portance of deep, high-resolution data to adequately characterize 
distant clusters. We were able to obtain a noticeable improve- 
ment on the previous X-ray characterization and our results set 
strong constraints on cluster evolution models. 

Our results are a first step towards understanding a possible 
cosmological evolution of the central feedback in galaxy clus- 
ters, and refining AGN feedback prescriptions in galaxy evolu- 
tion models at higher redshifts. To this aim, clearly a sizable rep- 
resentative sample of distant clusters with similar data quality, 
i.e. S/N and angular resolution, is needed. This is a task which 
only Chandra can currently achieve, albeit with a substantial in- 
vestment of time, whereas it is within comfortable reach of next 
generation wide field X-ray telescopes (e.g., the Wide Field X- 
ray Telescope, Giacconi et al. 120091 ). In Santos et al. ( 1201 II) . we 
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Studied how the measure of the parameter cs b varies with red- 
shift and angular resolution using simulated clusters. Based on 
this study we predict that we will be able to properly evaluate cs b 
for WARPJ1415 with an angular resolution of 5 "(Half Energy 
Width) or better. We conclude that only a combination of a large 
effective area and a large field of view, when coupled with a 
constant angular resolution of the order of 5", can provide a sig- 
nificant breakthrough on the study of cool-core clusters over a 
wide redshift range up to z=1.5, both in terms of statistics and 
data quality. With current X-ray facilities, we expect that only 
few cases can be studied with an accuracy comparable to that of 
WARPJ1415. 
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